Singularities in the optical response of cuprates. 
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We argue that the detailed analysis of the optical response in cuprate superconductors allows 
one to verify the magnetic scenario of superconductivity in cuprates, as for strong coupling charge 
carriers to antiferromagnetic spin fluctuations, the second derivative of optical conductivity should 
contain detectable singularities at 2A + A sp i n , 4A, and 2 A + 2A sp m, where A is the amplitude of 
the superconducting gap, and A s is the resonance energy of spin fluctuations measured in neutron 
scattering. We argue that there is a good chance that these singularities have already been detected 
in the experiments on optimally doped YBCO. 
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The pairing state in cuprate superconductors is pre- 
dominantly made out of Cooper pairs with d x 2_ y 2- 
symmetry [Q. This salient universal property of all 
high T c materials entails constraints on the microscopic 
mechanism of superconductivity. However, it does not 
uniquely determine it, leading to a quest for experiments 
which can identify " fingerprints" of a specific microscopic 
mechanism of d-wave superconductivity, a strategy sim- 
ilar to the one used in conventional superconductors |2j . 

Several resent experiments were interpreted as an indi- 
rect evidence that d x 2_ y 2 pairing in cuprates is produced 
by an exchange of collective spin fluctuations peaked at 
or near antiferromagnetic momentum Q = (n, 7r) In 
particular, the distance between the peak and the dip 
in the fermionic spectral function, A^oj), in ARPES ex- 
periments coincides with the frequency A s of the reso- 
nance peak measured in neutron scattering Q,[5| . This is 
exactly what one should expect for fermions interacting 
with a resonating spin collective mode [QJ^] (for phonon 
mediated superconductors, this is known as the Holstein 
effect ||). Similarly, a peak-dip structure of the SIS tun- 
neling conductance with peak-dip distance roughly con- 
sistent with A s has been obtained in the measurements 
on break junctions by Zasadzinski et al. for various dop- 
ing values (]]]. Carbotte et al. Q analyzed optical con- 
ductivity cr(w) in magnetically mediated d-wave super- 
conductors and argued that A s can be extracted from 
the measurements of the second derivative of a (to). 

In this paper we reexamine the behavior of the opti- 
cal conductivity in superconductors with quasiparticles 
strongly coupled to their own collective spin modes. Our 
results partly agree and partly disagree with those by 
Carbotte et al. || (see below). The key prediction of 
this paper, however, is novel: we argue that by measur- 
ing the frequency derivatives of the conductivity, one can 
not only verify the magnetic scenario, but, in principle, 
also independently determine both A s and A in the same 
experiment. 

Our argument goes as follows. For a superconductor, 
the real part of the conductivity, a\(uj), has a 8- func- 



tional piece due to the presence of the superconducting 
condensate. A nonzero o\ (uj) at a finite frequency is 
only possible if fermions have a finite lifetime. More pre- 
cisely, one of the two fermions exited in a process causing 
the AC conductivity should have a finite scattering rate, 
while another should be able to propagate, i.e., its energy 
should be larger than A. For clean, phonon-mediated su- 
perconductors, there are two sources for fermionic decay 
(see Fig. 1). One is a direct four-fermion interaction, 
which yields a threshold in the imaginary part of the 
self-energy, S" (w), at ui = 3 A - the minimal energy nec- 
essary to pull all three fermions in the final state out of 
the condensate of Cooper pairs. Another is the inter- 
action between an electron and an optical phonon. It 
yields the onset of S" (uj) at uj — A + £l p , where £l p is 
the frequency of an optical phonon Q (for simplicity, we 
assumed that the phonon propagator has a single pole). 
For the values of the coupling constant used to interpret 
the tunneling data in strongly coupled conventional su- 
perconductors like Pb ||, Q p > 2A, i.e., the onset of 
conductivity is at 3A + A = 4A (2A for dirty super- 
conductors Jl(| ) , while the signatures of phonon-assisted 
damping only show up at higher frequencies, and are un- 
correlated with the behavior of (Ti(uj) near 4A. 




FIG. 1. a) The exchange diagram for boson mediated in- 
teraction. The solid line stands for a propagating fermion. 
The wiggled line is a phonon propagator in case of electron- 
phonon interaction, and a magnon line in case of spin- fluctu- 
ation mediated interaction, b) The lowest order diagram for 
the fermionic self energy due to a direct four fermion interac- 
tion, also represented by a wiggly line 

For spin-mediated superconductivity, the situation is 
different. In the one-band model for cuprates, which we 
adopt, the underlying interaction is solely a Hubbard- 
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type four-fermion interaction. Spin excitations appear 
as collective modes of fermions, and their velocity v s is 
comparable to vp. For v s ~ vp, the low frequency spin 
dynamics is dominated by a decay process into a particle- 
hole pair and is purely relaxational in the normal state, 
with nearly featureless x"(Q,w) 0. (for phonon super- 
conductors the relaxation is also present but is strongly 
reduced due to a smallness of the sound velocity com- 
pared to vf p"l|]). 

Below T c , fermions acquire a gap, and spin-decay be- 
comes impossible for energies below 2A. The direct 
four-fermion interaction process in Fig. |l|b then yields 
a threshold in S" at 3A which gives rise to a singu- 
larity in the conductivity at w = 4 A jjjj. If x"(Q> w ) 
remained featureless, this would be the only effect. How- 
ever, several authors have demonstrated that the residual 
attraction in a d-wave superconductor binds a particle 
and a hole into a spin exciton at an energy A s < 2A. 
This effect gives rise to a peak in x"(Q,cl>) at u> = A s 
and makes it look like the spectral function for optical 
phonons. Accordingly, the conductivity acquires another 
threshold at 2 A + A s . Formally, this is analogous to the 
phonon case, but in distinction to phonons, A s < 2A. 
Then 2 A + A s < 4A, i.e., in clean systems, the lower 
threshold corresponds to the scattering by a spin exci- 
ton. Moreover, since both effects are due to the same 
underlying interaction, (the diagram in Fig. [l]b is just 
the first term in the series of graphs which constitute the 
spin-mediated scattering process shown in Fig [Ijb) the 
ratio A s /A and the relative intensity of the singularities 
in <Ti(u>) at 4A and 2A + A s are correlated. This correla- 
tion is a "fingerprint" of the spin-fluctuation mechanism. 
We will argue that there are strong indications that both 
singularities have been observed in the measurements of 
the optical conductivity in YBCO M , and their position 
and relative intensity are in reasonable agreement with 
the theory. 

Before we proceed with the calculations, a comment is 
in order. In the above discussion we neglected the mo- 
mentum dependence of the d-wave gap. Meanwhile, the 
computations of the optical conductivity involve averag- 
ing of the lifetime over the Fermi surface It is 
then a 'priori unclear whether the angular dependence of 
the d-wave gap with A(6) oc cos 2$ affects the positions 
of the two thresholds in the conductivity. Carbotte et al 
argued || that it does, and the singularity at 2A + A s 
(which they only considered) is determined by some av- 
eraged \2A(6)\ « A. We argue that averaging reduces 
strengths of the singularities but doesn't shift their po- 
sitions. Our argument is two-fold. First, we explicitly 
demonstrate below that the singularity in the conductiv- 
ity occurs at a frequency equal to the maximum value 
of the gap. Second, two of us and Finkel'stein argued 
earlier fLq ] that for spin-mediated d-w&ve superconduc- 
tivity, A(0) is at its maximum at hot spots (points at the 
Fermi surface separated by Q) JfE]. These are precisely 



the Fermi points which determine the position of the ex- 
citonic pole in x"(Q,w). Accordingly, the singularity in 
conductivity entirely comes from fermions near hot spots, 
and the threshold frequency 2A + A s involves a maxi- 
mum value of the gap and the resonance spin frequency 
at momentum Q. The same argumentation implies that 
4A threshold also involves a maximum value of the gap. 
Note for clarification that we are only considering here 
the singularities in the conductivity at u> > 2A. The reg- 
ular part of o~\(ijj) is not necessary confined to hot spots. 
In particular, for u> <C A the optical response is dom- 
inated by nodal quasiparticles for which E" is nonzero 
down to the lowest frequencies. 

We now proceed with the calculations. The real part 
of the optical conductivity in a superconductor is given 
by 



(Ti(uj) = Re- 



iS 



cMI CT (0,w) 



(1) 



where H a {9,uj) is the fully renormalized current-current 
correlator. In Matsubara frequencies, it is given by 
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and the normal and anomalous Green's functions are 



G k (iw m ) 
Fk(iu> m ) 



Sk(iWrn) + £k 



£ k (iw m ) - $l(iuj m ) - e k ' 

$k(^m) 

S k (iw TO ) - $ k (iw TO ) - e k ' 



(2) 

(3) 
(4) 



(we adsorbed a bare iui m term into £ k (ia-> m ). In prin- 
ciple, n CT is modified by vertex corrections related to 
dS/dfc jn]], but for spin- mediated scattering these cor- 
rections are small (see below). 

As an input for the computation of we need 
the forms of the fcrmionic self-energy E k (iw m ) and the 
anomalous vertex <& k (iw m ). We obtained these forms 
in Ref. |l6| by deriving and solving a set of Eliashberg 
equations within the spin-fermion model. This model 
adequately describes the interaction between low-energy 
fermions and their collective spin degrees of freedom 
fl P^|Jl6| at energies smaller than Ep- The full dynamical 
spin susceptibility peaked at (or near) Q mediates d x 2_ y 2 
superconductivity. As discussed, this susceptibility is by 
itself affected by low-energy fermions via a decay pro- 
cess into a particle and a hole, and has to be computed 
together with the fermionic self-energy and the pairing 
vertex. 

The justification of the Eliashberg approach for 
the spin-mediated superconductivity was discussed ear- 
lier |^,|l^,^6|, and we just quote the result: at strong di- 
mensionless spin-fermion coupling A, vertex corrections 
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and Vp dH / dk± , where k± is the component of the mo- 
mentum transverse to the Fermi surface, are small com- 
pared to dT./duj by log A/ A. In what follows we will ne- 
glect these corrections, i.e., approximate Ek(iw m ) and 
$k(«Wm) by E k (?w m ) = Y,{iui m ,9) and $(w m ,d). 
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FIG. 2. Frequency dependence of the real part of the 
optical conductivity <ri(w) at T = computed using the self 
energy and the pairing vertex determined from the Eliashberg 
equations for A = 1. The onset of the optical response is 
a) — 2 A + A s . The contributions from nodal regions (not 
included in calculations) yield a nonzero conductivity at all 
ui. They also soften the singularity at uj — 2 A + A s , but do 
not eliminate it. Insert - the behavior of the relaxation rate 
l/r(w) = (An/ujpi)Rell/cr(u))]. The frequency is measured in 
units Q which sets the overall energy scale in the Eliashberg 
solution. For A = 1, A = 0.204u> and A s = 0.291a>. 

As our goal is to study the singularities in ai(u>), we 
first perform calculations assuming that S and $ are in- 
dependent on 9 (i.e. that the superconducting gap is 
flat near the hot spots), and then analyze the results for 
a true cZ-wave gap. For a flat gap, the momentum in- 
tegration in Eq. (^) is straightforward. Substituting k 
integration by integration over and performing it, we 
obtain at T = and uj ^ 



n CT (iuj„) oc / duj' m d9 



D+D„(D + + D_) 



(5) 



Here E± = Z(iu>±,6), $± = $(iu)±,6), and D± = 
($| - S±) 1/2 , where w± = J ± lu/2. The conductiv- 
ity is obtained by converting this expression to the real 
axis ]T7|] . The singular piece in ui(w) near 2 A + A s can 
be obtained without a precise knowledge of £(o>) and 
the only information we need is that in a <i-wave 
superconductor, x"(Q,w) has a ^-functional singularity 
at a; = A s . This is what we found solving a set of three 
Eliashberg equations. Using this as an input and apply- 
ing a spectral representation for S" and <&" (which for a 
given x"{Qi w ) are described by a conventional set of two 
Eliashberg equations), we obtain that and $"(cj) 

are zero up to ui = A + A s , and undergo finite jumps 
at this frequency. By Kramers-Kronig relation, £' and 
<&' diverge as | log(w — where ujq = A + A s . The 



prefactor is the same for £' and Substituting these 
forms of S(a>) and $(w) into (0) we obtain after simple 
algebra that the conductivity emerges above 2A + A s as 
e 1//2 /log 2 e, where e = oj — (2A + A s ). This singularity 
obviously causes a divergence in the derivatives of the 
conductivity at e = +0. 

In Figj2| we show the result for the conductivity ob- 
tained by numerically solving Eq.(|^) using £(a>) and 
from Ref. [|l6|. We clearly see the expected thresh- 
old at 2 A + A s . The insert shows the behavior of the 
relaxation rate \/t{lo) — (An/Lu^Rell/a^)] where cj p i 
is the plasma frequency. Observe that \/t(uo) is linear in 
ui over a rather wide frequency range. This agrees with 
the earlier study of the normal state conductivity pq| . 
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FIG. 3. A calculated frequency dependence of 

W(lu) = ^[wRe[l/er(w)]] at T -> 0. This quantity is a sen- 
sitive measure of fine structures in the optical response. The 
locations of the extrema are: 1-2 A + A s , 2-4A, 3-2 A + 2A S . 
Dashed lines are the results are higher T. Observe that the 
maximum shifts to a lower temperature, but the minimum 
remains at 2A + A s . Inset- Experimental results for W(u>) 
at low T from Ref. The position of the deep minimum 
agrees well with 2A + A s . The extrema at higher frequen- 
cies are consistent with 4A and 2 (A + A 3 ) predicted by the 
theory. 

We next demonstrate that the position of the singular- 
ity is not affected by the angular dependence of the gap. 
Indeed, let the maximum value of the gap correspond to 
9 = and symmetry related points. At deviations from 
9 = both A and A s decrease. The decrease of A is ob- 
vious, the decrease of A s is due to the fact that resonance 
is a feedback from superconductivity, and its frequency 
scales as (A(9)) 1 /' 2 . Since both A and A s are maximal 
at a hot spot, we can expand u>o(0) = A(0) + A s (9) as 
ujq(9) — ujq — a9 2 , where a > 0. The singular pieces 
in E(w) and $(oj) then behave as | log(ojo — u> — a9 2 )\. 
Substituting these forms into @) and integrating over 9, 
we find that the conductivity itself and its first derivative 
are continuous at u = 2A + A s , but the second derivative 
of the conductivity diverges as d 2 a/duj 2 oc l/(|e| log 2 e) 
where, we remind, e = uj — (2 A + A s ). We see that the 
singularity is weakened by the angular dependence of the 
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gap, but it is still located at exactly 2 A + A s . 

The same reasoning is also applied to a region near 4 A. 
We found that the singularity at 4A is also weakened by 
the angular dependence of the gap, but is not shifted 
and still should show up in the second derivative of the 
conductivity. 

We now discuss the second derivative of the conduc- 
tivity in more detail. In Fig.|^ we present our numerical 
results for W{uo) = ^j(wRecr _1 (cj)) which effectively 
measures second derivative of conductivity (we followed 
Ref ||] and used the same W(lo) as for phonon supercon- 
ductors). We clearly see that there is a sharp maximum 
in W{uj) near 2A + A s followed by a deep minimum. We 
also see that W(ui) has extra extrema at 4A and, also, 
at 2loq = 2 A + 2A S . The last peak is a secondary effect 
due to a singularity in at u> — luq: <Ji(uj) is singular 
when the frequencies of both fermions in the polarization 
bubble exceed ujq. 

The experimental result for W(uj) in YBCO is shown 
in the insert. We see that the theoretical and experi- 
mental plots of W(oj) look rather similar, and the rel- 
ative intensities of the peaks arc at least qualitatively 
consistent with the theory. By the reasons which we dis- 
play below, we identify 2A + A s with the deep minimum 
in W(lo). This yields 2 A + A s w lOOmeV. Identifying 
the extra extrema in the experimental W(u>) with 4A 
and 2 A + 2A S , respectively, we obtain 4A ~ 130meV, 
and 2A + 2A S ~ 150meV. We see that three sets 
of data are self-consistent and yield A ~ 30meV and 
A s ~ 40 — 45meV. The value of A is in good agree- 
ment with tunneling measurements pofl , and A s agrees 
well with the resonance frequency extracted from neutron 
measurements pl| . We caution, however, that determi- 
nation of a second derivative of a measured quantity is 
a very subtle procedure. The good agreement between 
our theory and the experiment is promising but have to 
be verified in further experimental studies. Nevertheless, 
our calculation clearly demonstrates the presence and ob- 
servability of these " higher harmonics" of the optical re- 
sponse at 4 A and 2 A + 2A S . 

So far we considered only the singular part of o~\(uj). 
In Fig.^j we compare our results for o-\{uj) (ignoring the 
contributions from the nodes) directly with the experi- 
mental data by Puchkov et al. for optimally doped 
YBa 2 Cu 3 06+5. We used A = 1 and the overall energy 
scale uj which yield A ~ 30meV and A s ~ 45meV 
as the solution of the Eliashberg set, and also oj p = 
1.2 x 10 4 cm _1 , similar to that in |lS| ], We see that the fre- 
quency dependence of the conductivity at high frequen- 
cies agrees well with the data. The measured conduc- 
tivity drops at about lOOmeV in rough agreement with 
2A + A s « lOOmeV in our theory. As in earlier stud- 
ies [ |b4||i"g| ] , to match the magnitude of the conductivity, 
we had to add the constant 7 x 10 _4 O cm to (<7i(a;)) _1 . 
We view the good agreement between theory and experi- 



ment at lu > 2A + A s is predominantly an indication that 
the momentum dependence of the fermionic dynamics be- 
comes irrelevant at high frequencies, and fermions from 
all over the Fermi surface behave as if they were at hot 
spots. The insert to Fig.|| shows cr^ 1 (ix)). We see that it 
is linear over a substantial frequency range. 
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FIG. 4. A comparison of the theoretical result with the 
experimental data of Puchkov et al. Iljj. The substructure 
in the theoretical o\ at low frequencies is an artifact. Insert - 
the behavior of a^ l {ui). 

We emphasize, however, that this linearity is only an 
intermediate asymptotic. At the highest u>, our theory 
yields cri(w) oc oj -1 / 2 . The lower boundary for a; -1 / 2 be- 
havior decreases with increasing A. For optimally doped 
Bi2212, A is somewhat larger than in YBCO as A « A s , 
and we expect a more pronounced uj^ 1 / 2 behavior at high 
frequencies. This trend is consistent with the data of 
Ref. P^| . This issue, however, requires further study as 
<7\{uj) cc uj~ x / 2 at intermediate frequencies was also ob- 
tained in Ref. [jl4| assuming a strongly momentum de- 
pendent scattering rate. 

Finally, we comment on the position of the 2A + A s 
peak and compare our results with those by Carbotte 
et al. Theoretically, at T = and in clean limit, 
the maximum and minimum in W(lo) are at the same 
frequency. We found, however, that at finite T, they 
quickly move apart. We present the theoretical temper- 
ature dependence of W{lu) in Fig ||[ Carbotte et al. || 
focused on the maximum in W{uj) and argued that it is 
located at A + A s instead of 2A + A s . We also found that 
the maximum in W(uj) shifts to a lower frequency with 
increasing temperature, already at T where the temper- 
ature dependence of the gap may be neglected. On the 
other hand, the minimum in W(u>) moves very little with 
increasing T and virtually remains at the same frequency 
as at T = 0. This is our reasoning to use the minimum in 
W{ui) as a much more reliable feature for the comparison 
with experiments. This reasoning is in agreement with 
recent conductivity data on optimally doped Bi2212 |2^] 
- W(lo) extracted from these data shows strong down- 
turn variation of the maximum in W(lo) with increasing 
temperature, but the minimum in W(uj) is located at 
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around HOmeU for all temperatures. 

Finally, we briefly consider whether one can extract a 
resonance spin frequency from the measurements of the 
Raman intensity in a d— wave superconductor. The Ra- 
man intensity is given by p3[ 



R(ui) =Im J d6V 2 {6)n R (u,( 



(6) 



where V{&) is Raman matrix element, and is same 
bubble as for conductivity, but with a different sign of the 
anomalous FF term. The latter is a consequence of the 
fact that Raman vertices are scalar and do not change 
sign under k— > —k. Performing the integration over 
quasiparticle energies in the same way as for conductivity 
we obtain in Matsubara frequencies ]2q| 



Ii a (iuj n ) oc / duj' m dO- 



D+D- 



D+D_(D + + D_) 



(7) 



For mostly studied B± g scattering, the Raman vertex has 
the same angular dependence as the <i-wave gap, i.e., 
V(9) oc cos (29) [ p3p^ ]. Straightforward computations 
then show that for a d— wave gas, R(u) oc ui 3 at low 
frequencies [M, and diverges as to approaches 2 A first 
as \u> — 2A| 1 '^, and then, in the immediate vicinity of 
2A, as log|w - 2A| At larger frequencies R(u) 

gradually decreases. At strong coupling, we performed 
the same analysis as for conductivity and found that the 
sign change of the FF term in the bubble, compared to 
that for conductivity, has a drastic effect: near 2A + A s , 
singular contributions from E + I]_ and <& + $_ terms in 
Eq.(^) cancel each other. As a result, for a flat gap, only 
the second derivative of R(oj) diverges at 2A + A s . For a 
quadratic variation of a gap near its maximum, the sin- 
gularity is even weaker and shows up only in the third 
derivative of R(uj). Obviously, this is a much weaker 
effect than that for conductivity, and its determination 
requires a high quality of the experiment. Notice, how- 
ever, that due to the closeness of hot spots to (0, 7r) and 
related points, at which vf vanishes, the actual smearing 
of the singularity due to the angular integration may be 
less drastic than in our theory and the singularity in R(uj) 
at 2A + A s may possibly be extracted from the data. 

To conclude, in this paper we examined the singular- 
ities in the optical conductivity in a d-wave supercon- 
ductors assuming that the pairing is mediated by over- 
damped spin fluctuations. We argued that a\(uj) should 
have singularities at 2A + A s , 4 A and 2A + 2A S , where A 
is the maximum value of the d— wave gap, and A s < 2A 
is the resonance spin frequency. The experimental detec- 
tion of these singularities would be a strong argument in 
favor of the magnetic scenario. We argued that there is a 
good possibility that all three singularities have actually 
been detected in recent data on YBCO. 
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